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Josephson tunnel junctions are proposed as efficient phase-tunable thermal rectifiers. The latter exploit the strong tem- 
perature dependence of the superconducting density of states and phase-dependence of heat currents flowing through 
Josephson junctions to operate. Remarkably, large heat rectification coefficients up to ~ 800% can potentially be 
achieved using conventional materials and standard fabrication methods. In addition, these devices allow for the in-situ 
fine tuning of the thermal rectification magnitude and direction. 
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Electronic circuits consist of a number of components (e.g., 
transistors, diodes and switches) connected together to enable 
the execution of different operations. In metals, electrons are 
responsible of energy transport as well, in what is commonly 
referred to as electronic heat transport. ^'^ It is therefore natural 
to address the feasibility of networks that might eventually al- 
low for the implementation of, for instance, thermal computa- 
tion, thermal logic operations or data storage (for a review see 
Ref. 3 and references therein). Additionally, mastering of heat 
currents represents an important breakthrough in different re- 
search fields of nanoscience such as solid-state cooling,^ ra- 
diation detection ^ quantum computing^ or the emerging field 
of coherent caloritronicsr'~^ So far, a strong effort has been 
devoted to envision thermal rectifiers, i.e., structures allow- 
ing high heat conduction along one direction but suppressed 
thermal transport upon temperature bias reversal. ^'^ Most of 
these proposals deal with phononic heat transport, ^^~^^ very 
few deal with electronic heat conduction ^^~^^ and even less 
have demonstrated feasible experimental realizations. ^^'^^ 

In this Letter we propose and analyze theoretically the per- 
formance of a thermal diode consisting of a SIS' Josephson 
tunnel junction, where the I stands for an insulating barrier 
and S and S' represent two different superconducting elec- 
trodes. Although never considered so far for such a pur- 
pose, superconducting tunnel junctions appear particularly 
well suited for the implementation of electron heat rectifiers. 
Heat transport in such structures is deeply influenced by the 
strong temperature dependence of the superconducting den- 
sity of states (DOS). Yet, the Josephson effect provides the 
thermal diode with an even more interesting capability. As re- 
cently demonstrated,^ heat currents flowing through Joseph- 
son tunnel junctions depend also on the macroscopic quantum 
phase difference of the Cooper pair condensates, ^^~^^ just as 
charge currents do. This latter property has outstanding con- 
sequences enabling to conceive thermal quantum devices that 
go well beyond the simple concept of heat rectification. All 
these features lead, under suitable conditions, to remarkable 
rectification coefficients as large as ~ 800%. 

We shall start, first of all, by defining a heat rectification pa- 
rameter IR. To this end, let us consider two different supercon- 
ductors. Si and S2, weakly coupled so to implement a Joseph- 
son tunnel junction as shown in Fig. 1(a) and 1(b). Each 
superconductor is characterized by its energy gap Ai and A2 




FIG. 1. (a) and (b) Josephson thermal diode scheme correspond- 
ing to the forward and reverse thermal bias configuration, respec- 
tively. The diode is phase biased (9) and temperature biased with 
Ti ^ T2 as well, but the voltage across the junction vanishes, (c) 
'51 vs. Thot for different values of r^oid corresponding to a normal 
metal/insulator/superconductor (NIS) thermal diode so that A2 = 0. 



leading to critical temperatures T^i and 1^2, respectively. The 
electronic temperature in both Si and S2 is kept at fixed T\ 
and 72, respectively, and the voltage drop across the junction 
is set to zero. Additionally, (p denotes the macroscopic phase 
difference across the junction with normal-state resistance 7?j. 
In the forward thermal bias configuration, a thermal gradient 
is created by setting T\ = T^ot > T2 = T^oid which leads to a 
total heat current /+ flowing from Si to S2 [see Fig. 1(a)]. In 
the reverse thermal bias configuration, the thermal gradient is 
inverted so that 7i = Tcoid < T2 = Thot leading to a heat cur- 
rent /_ flowing from S2 to Si [see Fig. 1(b)]. Under these 
hypothesis we define the rectification coefficient as 



:r(%) : 



/_ 



xlOO. 



(1) 



We describe now the equations governing heat transport in 
the Josephson thermal rectifier. We focus on the electronic 
contribution to heat transport only, and neglect eventual heat 
currents carried by lattice phonons. The forward and reverse 
total heat currents flowing through the Josephson junction 
read^ 



J^-) ='^^p[^hot(cold),^cold(hot)] 
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-Jint [^hot(cold), 7iold(hot)] COS (jf), 
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where the term Jqp accounts for the energy carried by 
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FIG. 2. Panels (a), (b) and (c) show three density plots of Jl vs. 
Thot and r calculated for cp = 0, (p = n/2 and (p = n, respectively. 
Panels (d), (e) and (f) show three selected profiles of Jl\s r for the 
same values of (p corresponding to the colored dotted lines in (a), (b) 
and (c). Notice that the scale is logarithmic above the break in the 
vertical axis. In addition, a dashed line indicates 31 = 0. All curves 
have been calculated for Tcoid = O.OlTci. 



quasiparticles, 
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and k, I =hot,cold. 
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Eq. (3), J<a{£jk) = 

is the normalized smeared (by non 



V[e/A«(r^)+/7]2- 

zero 7) BCS quasiparticle DOS in S^ with a = 1,2.^"^ In the 
following we will assume 7~ 10~^, which describes realistic 
SIS' junctions.2^-2^ Additionally, f{Tk) = (1 + e^/^^^k)'^ 
is the Fermi energy distribution, Aa (Tj^) is the temperature- 
dependent energy gap of S^, ^b is the Boltzmann constant 
and e is the electron charge. On the other hand, the term Jtnt 
in Eq. (2) is the amplitude of the phase-dependent component 
of the heat current,^^"^^ 



e^Rj Jo 



(4) 



mThot)-f{Tcold)]d£, 



whereMa(e,7i)=3 
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This component 



is peculiar to the Josephson effect and arises as a consequence 
of tunneling processes through the junction involving both 
quasiparticles and Cooper pairs. ^^~^^ It is worthwhile to em- 
phasize that, depending on the value of (p, the second term in 
Eq. (2) may switch its sign. The existence and sign of this 
component was experimentally demonstrated in Ref . 5 

It is illustrative to start by analyzing the case in which one 
of the two electrodes is a normal metal, i.e., a NIS junction. 



FIG. 3. (a) and (b) Jl vs. Thot for a few values of Tcoid, at (jp = and 
(p = 7i, respectively. Notice that in (b) the vertical scale is logarithmic 
for Jl> 0. (c) Periodicity of 31 with (p calculated for the same values 
of Tcoid and the corresponding Thot that maximize Jl. The dashed line 
indicates 31 = 0. All these curves correspond to r = 0.75. 



For this purpose we can simply set A2 = which leads to the 
complete suppression of Jtnt in Eq. (2). We calculate Jl as 
a function of Thot for different values of Tcoid, i.e., for dif- 
ferent temperature gradients established across the weak link. 
As shown in Fig. 1(c) a maximum positive rectification of 
- 26% is obtained for Thot - O.SSJd at T^oid = O.Oird. As 
Thot increases, heat rectification starts to decrease eventually 
inverting its sign which implies that heat flux from the normal 
metal to Si becomes preferred. Furthermore, by increasing 
Tcoid leads to a reduction of 31 which reaches its maximum for 
larger values of Thot- We note than, a bare NIS junction offers, 
in this way, an already quite rich response in terms of heat 
rectification. 

In the following analysis we shall focus on the case in 
which both junction electrodes are superconductors, i.e., a 
SIS' junction. We define for clarity r = A1/A2 < 1. By fix- 
ing the temperature of the second electrode to Tcoid = O.OlTd, 
we calculate IR as a function of Thot and as a function of r. 
The result is plotted in Fig. 2(a), (b) and (c) for three repre- 
sentative cases, corresponding to 9 = 0, 9 = 7r/2 and (p = 7i, 
respectively. Three selected profiles of ^ as a function of r 
for different Thot values are shown as well in Fig. 2(d), (e) and 
(f) for the same values of (p. 3i depends strongly on r reach- 
ing its maximum at r c^ 0.75 and Thot — O.llTd for (p = 71, 
dropping then to zero at r = 1 . The inspection of these graphs 
also reveals how phase biasing across the junction does make 
a substantial difference. In particular, the heat rectification co- 
efficient does not only change by almost two orders of magni- 
tude from (p = 0to (p = 71 but it also switches its sign. 

The dependence of 31 on Thot for (p = is shown in Fig. 
3(a) for different values of Tcoid and r = 0.75. In particular, 
31 is negative as long as Thot < T^, i.e., as long as one of the 
two electrodes remains in the superconducting state, reaching 



values ofJlc^ -25% at Thot ^ 0.767^1 and T^oid = 0.017^1. 
Above this point, the sign of Jl and, therefore, the rectifica- 
tion direction depends on the thermal gradient being however 
slightly reduced. On the other hand, curves corresponding to 
9 = TT are shown in Fig. 3(b) vs. T^ot for the same values of 
Tcoid. Remarkably, large values of IR ~ 800% can be obtained 
at Thot — O.llTci and Tcoid = O.OlTd. If Tcoid is increased, ^ 
is in general reduced reaching its maximum for larger values 
of Thot- Above T^, the behavior is identical to that calculated 
for 9 = since one electrode remains always in the normal 
state and the phase plays no role any more. In Fig. 3(c), the 
periodicity of Jl with (p is shown for the same values of Tcoid 
and the corresponding Thot that maximize Jl. We stress how 
(p-dependence provides the Josephson thermal diode with an 
unique tunability. 

It is worthwhile to emphasize that the SIS' junction rectifies 
heat only if Ai ^ A2. As for the case of the NIS diode, heat 
rectification demands the combination of two different DOS 
being (at least one of them) strongly temperature-dependent.^^ 
Yet, a fairly large thermal gradient is required as well. Indeed, 
heat rectification is absent in the linear-response regime, that 
is, for small temperature differences 5T = Thot — ^coid «T = 
^°^^ ^°^^ . In this case, the total electron heat current flowing 
through the Josephson junction reduces to^^ 
5T 
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'le^hT^RjJo \2kBTJ (5) 

X [1^1 (e, r)?^2(e, T) - Ml (e, r)M2(e, T) cos (p] de, 

which depends on the average temperature of the two elec- 
trodes only. 

We discuss in the following some possible experimental re- 
alizations faced to make the most of the phase-dependence 
of the Josephson thermal rectifier. On the very first place we 
are interested in playing with the macroscopic phase differ- 
ence across the Josephson junction during operation. Phase 
biasing of a Josephson junction can be achieved, in gen- 
eral, through supercurrent injection or by applying an exter- 
nal magnetic flux. In the former case, schematized in Fig. 
4(a), a Josephson current /j is forced to flow through the junc- 
tion via two extra control superconducting wires connected to 
the diode's core through clean or tunnel contacts. The con- 
trol wires can be made of a third superconductor S3 with 
energy gap A3 >> Ai,A2 so to suppress heat losses. The 
phase-current relation under such circumstances is given by 
sin (9) = i]/iy^^ Provided that /j < /j, /j being the junction 
critical current, a phase gradient contained within the sec- 
tions —7i/2 < (p < 71/ 2 can be established. An analogous 
phase gradient can be obtained using a direct current super- 
conducting quantum interference device (DC SQUID) pierced 
by an external control magnetic flux <I> as shown in Fig. 
4(b). If both junctions are identical and neglecting the loop's 
geometrical inductance, the phase-flux relatio n must satisfy 
cos{(pa) = cos{(piy) = Y^^Il~^^cos(2^r^7^o)I7^ where % and 
(ph are the phase drops across each junction^'^ and <I>o is the 
flux quantum. The optimum phase configuration in terms of 
heat rectification, i.e., (p = n, can be reached by using a rf 
SQUID as shown in Fig. 4(c). Fur such a purpose, the ther- 
mal diode can be enclosed through clean contacts within a 




FIG. 4. Panels (a), (b) and (c) show three possible ways of in-situ 
phase biasing of a Josephson junction within —71/ 2 < (p < n/2 [(a) 
and (b)] or within —K<(p<K (c). (d) Possible experimental real- 
ization of the Josephson thermal rectifier. Source and drain normal 
metal electrodes are tunnel-coupled to the diode's core. Additional 
superconducting probes tunnel-coupled to source and drain allow for 
the implementation of SINIS thermometers and heaters. 



superconducting ring S3 pierced by a control flux <I>. Ne- 
glecting again the loop's inductance, the phase-flux relation 
is given in this case by 9 = 2tc^/^q^^ enabling the phase 
drop across the junction to vary within the whole phase space, 
i.e., -n <(p <n. 

Figure 4(d) shows a device envisioned to probe experimen- 
tally the effects discussed above. Two identical normal metal 
electrodes, source and drain, are weakly connected one each 
via a resistance R^ to both Si and S2, respectively. Supercon- 
ducting probes can be tunnel-coupled to these electrodes so to 
implement SINIS thermometers and heaters.^ Yet, the forward 
thermal bias configuration can be realized by intentionally in- 
creasing the electronic temperature in source electrode up to 
^src = ^h and probing the temperature in drain electrode T^ . 
On the reverse configuration, we set T^ = Th and T~^ is mea- 



sured in a similar way. The difference 57 e = T^ — T~^ for 
a given Th can be used to assess experimentally heat rectifi- 
cation. 5% can be computed numerically by solving a sys- 
tem of thermal equations accounting for the heat exchange 
mechanisms present in our device [see Fig. 5(a)]. On the for- 
ward configuration, electrons in Si exchange heat with elec- 
trons in the source at power /src^Si and, at power /+, with 
electrons in S2. On the other hand, electrons in S2 exchange 
heat with drain's electrons at power /s2^dr- Finally, electrons 
in the whole structure exchange heat at power /g-ph with lat- 
tice phonons that we assume to reside at bath temperature 
^bath- Under such circumstances, the three unknown quanti- 
ties, i.e., Thot, ^coid and T^ can be calculated for given initial 
conditions by solving the following system of thermal balance 
equations 
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-^src^Si (^h, ^hot) - '^e-ph,Si (^hot) - J^ (^hot, ^cold) = 0, 
J^ {Thot, ^cold) - '4-ph,S2 (^cold) - Js2^dY{Tco\d, T^) =^i (6) 
-^82^(11(^00^ 7 ^dr ) ~ -^e-ph^rl^bath, ^dr ) ~ ^- 
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FIG. 5. (a) Thermal model describing the main heat exchange mech- 
anisms playing a role in the device proposed in Fig. 4(d). Panels 
(b) and (c) show the computed values of IR as a function of T\^ and 
(p, respectively. Dashed lines in (b) correspond to the curves in (c). 
Panels (d) and (e) show the same curves corresponding to S7e- 



In the above expressions, /e-ph,src(dr) = ^N^N{T^rcidv) " ^blth)^' 
Vn and En being the volume of the normal metal electrode 
and the electron-phonon coupling constant, respectively. Fur- 
thermore, we assume Tbath = 10 mK << Tj^ « A(r^)/^B 

SO that /e-ph,s« - 0.95i:sVs7;^^ Vk ,3i with Vs and Zs be- 
ing the volume of each superconducting electrode and the 
electron-phonon coupling constant, respectively. As represen- 
tative parameters we set Vn = Vs = 2 x 10~^^ m^, R] = 10 
k^ and R^ = 100 Q.. Source and drain electrodes can be 
made, for instance, of Cu for which Zn — 3 x 10^ WK~^m~^^ 
whereas the diode's core can be made of Al and Mn-doped Al 
with Ls — 0.3 X 10^ WK~^m~^,^ since the latter allows for 
fine tuning of the aluminum superconducting gap.^^'^^ In this 
= 1.4 K and r = 0.75. With this 



way, we set T^i = jj^^ 
set of parameters, ^ can be determined as a function of Th 
for given values of (p. The resulting curves are plotted in Fig. 
4(b) together with the computed values of 57e vs. 7h which 
are plotted in Fig. 4(d). Remarkably, in the present setup, a 
maximum Ji ~ 340% can be reached. The latter corresponds 
to a temperature difference as large as 5 T^ ~ 140 mK which 
is easily measurable with standard SINIS or SNS thermome- 
try techniques. ^'^ Even more interesting, phase-coherence fin- 
gerprints are clearly observable as well. Notably, Ji and 5% 
show the expected 2;r-periodicity as shown in Fig. 4(c) and 
(e). In closing, we emphasize that, at such low bath temper- 
atures, both superconducting electrodes are only marginally 
coupled to phonons.^ Indeed, neglecting the contribution of 
'^e-ph,s, leads to differences less than ~ 5% of the values pre- 
sented here. 



To summarize, we have proposed and analyzed the concept 
of a Josephson thermal rectifier. Under appropriate condi- 
tions, a remarkably large rectification coefficient of ^ ~ 800% 
can be obtained. In addition, the Josephson thermal diode 
is phase- tunable. This latter property allows to maximize 
heat rectification in-situ or, even, to switch its sign. Such 
a device might find a straightforward application, e.g., in the 
field of electronic refrigeration enabling magnetic-flux depen- 
dent heat management and thermal isolation at the nanoscale. 
The operation principle which is at the basis of this heat rec- 
tifier will likely contribute, on the other hand, to improve 
the performance of other different coherent thermal compo- 
nents such as heat transistors or splitters. ^~^ These thermal 
devices might potentially lead to the emergence of coherent 
caloritronic nanocircuits. 
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